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We study theoretically the ac Josephson effect in voltage biased planar junctions 
of d-wave superconductors. For some orientations of the superconductors a current 
peak is found at finite voltage in the current-voltage characteristics. We pick out 
the relevant physical processes and write down an analytical formula for the current 
which clearly shows how the midgap state acts as a resonance and produces the 
peak. We present a possible explanation for the zero-bias conductance peak, recently 
found in experiments on grain boundary junctions of high-temperature supercon- 
ductors, in terms of resonant transmission through midgap state of quasiparticlcs 
undergoing multiple Andreev reflections. We note that within our framework the 
zero-bias conductance peak appears in rather transparent Josephson junctions of 
d-wave superconductors. 
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1. Introduction 

The controversy regarding the symmetry of the order parameter in high-temperature superconduc- 
tors (HTS) seems today, to some extent, have been resolved in favour of the d 2 .2_ a 2-wave (d-wave) 
symmetry jl], ^, Since current transport through superconducting junctions is sensitive to both 
magnitude and phase of the order parameter, experiments on Josephson junctions can measure the 
intrinsic phase of the d-wave order parameter. Therefore, a great deal of attention has been directed 
towards the Josephson effect in HTS junctions. Both the dc @, |, §, § and ac §, % |k| [l| 
effects have been studied. 

An important discovery was that at surfaces and interfaces of d-wave superconductors a midgap 
state (MGS) may be formed [ fill) , which affects the current transport. A quasiparticle reflected at 
the interface will sense different order parameters before and after the scattering event since the 
momentum is changed. When there is a sign difference between the order parameters before and 
after scattering, a zero-energy bound state is formed (MGS). In normal metal-d-wave superconduc- 
tor (N|d) junctions the MGS gives rise to a zero-bias conductance peak (ZBCP), established both 
theoretically |D| [y| and experimentally @ [H], |l9| . Previously, it has been shown that MGS 
gives rise to a current peak at finite voltage in both s-wave superconductor-d-wave superconduc- 
tor (s\d) and d-wave superconductor-d-wave superconductor (d|d) junctions, resulting in negative 
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Fig. 1. The model system we are considering. The orientations of the two ci-wave superconductors are given by the 
angles ol\ and c*2- An electron-like quasiparticle is incident on the junction at the angle 9. We have here introduced 
short normal regions on each side of the barrier. The strength of the barrier can be tuned from the ballistic to the 
insulating case. 



differential conductance §, |(J |l|, The main purpose of this paper is to analyze the results of 
Refs. H and jlO| in terms of the physical processes giving rise to the peak at finite voltage. 

Knowing the effects of MGS in the N\d junction, it is not straightforward to predict the effects 
in the s\d and d\d junctions. This is the case because the mechanisms involved in the current 
transport are more complicated when both electrodes are superconducting. In a voltage biased 
junction between two superconductors, the phase difference over the junction is time-dependent 
according to the Josephson relation <p = 2eV/h. This makes the scattering in the junction inelastic: 
a quasiparticle incident on the junction at energy E produces a scattering state with sideband 
energies E n = E + neV , where V is the voltage over the junction and n is an integer. 

The knowledge about current transport through voltage biased s-wave s\s junctions has greatly 
increased in recent years. Using a Landauer-Biittiker scattering method (extended to include su- 
perconducting electrodes) it has become possible to understand s\s junctions with arbitrary trans- 
parency of the insulator in terms of multiple Andreev reflections (2^, |2^, p5[ ^6). The theory 
has also been extended to include d-wave superconductors [|[ [h]] . This theory is used in the present 
paper in order to study the effects of MGS in voltage biased superconducting junctions. 

The rest of the paper is organized in the following way. In Chapter 2 we briefly describe the model 
and how to calculate the current-voltage (IV) characteristics. In Chapter 3, the scattering problem 
is solved analytically taking into account only processes giving rise to the current peak. We study 
the s\d junction (where the peak is most pronounced) and compare the current contribution from 
the processes responsible for the current peak with the complete (numerically calculated) IV-curve. 
In Chapter 4 we present a qualitative discussion of how a ZBCP may appear in voltage biased d\d 
junctions. This study was motivated by recent experiments |2^, ^J. In Chapter 5 we summarize the 
paper. 



2. The model 



The system we are considering is a planar junction between d-wave superconductors as shown in 
Fig. pi. The superconducting gap function depends on the quasiparticle's direction of propagation 9. 
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We assume that the superconductors have pure <i-wave symmetry and let A(0) = Ao cos[2(# — a)], 
where a is the orientation of the superconductor relative to the interface normal. For an s-wave 
superconductor, A(6>) = A s . 

We neglect surface roughness in our treatment of the junction region, modelling the barrier by 
the potential V(x) — HS(x). The reflection and transmission amplitudes of the barrier are therefore 
angle-dependent: r{9) = Z/(icos9 — Z) and t(0) = icos9/(icos9 — Z), where the dimensionless 
parameter Z — 2mH/h 2 describes the transparency of the junction ^7|. A small value of Z corre- 
sponds to a transparent junction and a large value corresponds to a tunnel junction. The limiting 
cases are Z — (the ballistic junction - no scattering) and Z — oo (two uncoupled electrodes). 
Having only specular reflections implies that the quasiparticle's momentum along the junction k y 
is a conserved quantity. This means that the current is an average over all injection angles where 
the current contribution from each injection angle can be calculated separately. The difference from 
a junction between s-wave superconductors is that we here must take into account the angle de- 
pendence of the (i-wave gap. In particular, we must remember that the gap changes after normal 
reflection at the barrier. 

Self-consistency of the gap is important in junctions of (i-superconductors p| . It has been shown 
in the tunnelling limit that including self-consistency of the gap produces bound states with non- 
zero energy along with the MGS. The length of the normal regions introduced on each side of the 
barrier in Fig. [l] could be thought of as a model of the possible suppression of the d-wave gap near 
the interface. Since we in this paper focus on effects of MGS, which is due to the intrinsic sign of 
the d-wave order parameter only, we neglect these effects and put the length of the normal regions 
to zero. We therefore assume a step-like behaviour of the gap-functions: 

Ai(0), x<0 

A 2 (6»)e^°, x>0. [) 

Since the overall phase is unimportant we can here choose Ai real and let the phase of the right-hand 
superconductor be equal to the phase-difference (f>o over the junction. 

The time-dependent Bogoliubov-de Gennes equation describing anisotropic superconductors is 
solved for a voltage biased junction, modeled as described above, by the method described in detail 
in Ref . jl(| . We obtain the scattering states by matching ansatz wavefunctions at the NS interfaces 
and at the barrier for all injection angles 9 and energies E. The current is a sum over contributions 
from all scattering states, where each contribution is found by inserting the wavefunction into the 
current formula. The sum is in the end turned into an integration over injection angle and energy. 

We choose to calculate the current in the normal region to the left of the barrier. The wave 
function in this region is 

where the side-band energy is E n = E + neV. The index a = {e^, h~ *, h^, e*~ } labels the four types 
of incoming quasiparticles: electron-like and hole-like quasiparticles injected from the left and right 
electrodes. Here is \i e / h the wave vector of electrons and holes respectively. In our two-dimensional 
problem we have k = (k x ,k y ) = k(cos9,smd) and k = (—k x ,k y ) — k(cos9,sm9), where 9 = it — 9 
is the angle after normal reflection. The wave function coefficients a, 6, c, and d in Eq. (^|) are found 
by solving the matching equations. 

The current per a6-plane for a junction of width L y is 

I 1 f 7 */ 2 f°° ^ 

- = 7m dOcosO dEf(E)J2j°(e,E), (3) 
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where 00 = 2ekpL y D /irh is the normal state conductance, D = J ddT>cos0/2 is the normal state 
transparency averaged over all angles 8 (V — \t(9)\ 2 ), Ep is the Fermi energy, and f(E) = 1/[1 + 
exp(_E/fceT)] is the Fermi distribution function. In Eq. we sum over electron-like and hole-like 
quasiparticles injected from the left and right reservoirs described by the index a. 

The current density appearing in Eq. (^) is expressed in terms of the wavefunction coefficients in 



f(9,E) = N°(6,E) £ K| 2 - K| 2 + K| 2 - K| 2 . (4) 



In order to collect contributions with the same momentum along the junction k y , we inject the 
four quasiparticles a at four different angles with the following four gaps: A a = {Ai, Ai, A 2 , A 2 } 
suppressing the dependence on 9. A consequence of this is that four density of states N"(9,E) = 
Q(\E\ - \A a (9)\)\E\/^E 2 - A°(8) 2 appear in Eq. (§. 

3. Midgap state resonance 

Consider an N|d^/4 junction biased at arbitrary small voltage V and an electron incident on 
the junction from the normal metal side at the Fermi energy. Due to MGS, the electron will be 
Andreev reflected with unit probability, independent of the strength of the interface barrier. In 
energy space, this can be thought of as resonant transmission through MGS: the electron at energy 
E is turned into a hole with energy E + 2eV with unit probability. In the Andreev reflection process 
a charge 2e is transferred into the superconductor (forming a Cooper pair) and a current flows 
through the junction. This current flow through MGS is the explanation of the ZBCP seen in many 
experiments 0, |T| |l9|. 



The advantage of the above description of current flow through MGS is that it can be gener- 
alized to the case when the normal metal is changed to a superconductor. In a junction between 
two superconductors an incident quasiparticle will create both electrons and holes in the normal 
region, which will undergo multiple Andreev reflections. Because of the voltage drop over the nor- 
mal region each particle changes its energy by eV each time it passes the normal region. In this 
way a scattering state with amplitudes at all sideband energies E n is built up |2^| . If we consider 
a quasiparticle incident from the left superconductor at energy E, Andreev reflections at the right 
superconductor will take place at the odd sideband energies E 2n +i, while Andreev reflections at the 
left superconductor will take place at the even sideband energies E 2n . When e.g. the right electrode 
is a d-wave superconductor, oriented in such a way that MGS is formed, a transmission resonance 
is produced in energy space: the probability of transmission from energy E 2n to energy E 2n+2 is 
unity if E 2n+1 = 0. 

The resonant transport through MGS described above results in current peaks in the IV- 
characteristics of junctions containing d-wave superconductors, as previously reported in Refs. || 
and JTo| . Since the current peak is most pronounced in the s|d,r/4 junction we will here concentrate 
on this particular junction. 

3.1. Midgap state resonance in the junction 

Generally, it is possible to express the dc current in superconducting junctions as a sum over 
contributions from different rt-particle processes, where n is an integer. In a single scattering state, 
originating from a quasiparticle incoming at energy E, the n-particle contribution to the dc current 
is n ■ 1%, where 1% is the outgoing probability current at energy E n = E + neV (^6). In this section 
we will focus on the resonant two-particle process involving the MGS and derive the expression for 
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Fig. 2. The two-particle process in the s|d w /4 junction involving MGS for a quasiparticle injected from the left 
electrode. When matching, we divide the scattering state into three parts: the part in between the injection point 
and the exit point (region I), the part below the injection point (region II), and the part above the exit point (region 
III). The MGS resonance is marked by the thick lines. 



its contribution to the current explicitly, without using the concept of probability current. To be 
specific we will calculate the dc current contribution from quasiparticles incoming from the left, at 
energy E ps — eV, that are Andreev reflected through the MGS on the right hand side, at E\ f» 
and then leave the normal region to the left at energy E 2 ~ eV (see Fig. ||). This two-particle 
process involves one Andreev reflection, implying a net current transport of 2e, as we will also see 
in the final expression for the current. 

In order to solve the scattering problem analytically, we divide the energy axis into three parts: 
the part in between the injection point and the exit point, the part below the injection point, and 
finally the part above the exit point. We denote these three parts by I, II, and III, as shown in Fig. [| 
The problem of calculating the scattering state can be mapped onto the problem of calculating the 
wavefunction for tunneling through a one-dimensional multi-barrier structure (in energy space), as 
drawn in Fig. ||. The structure of the scattering state is better illustrated by Fig. ||, which also makes 
clear the way the matching of the ansatz wave functions is done and how resonances may appear. 
From now on, the discussion will therefore be connected to Fig. [|. 

The coefficients in part I are coupled by a scattering matrix S 2 o'- 



"2 



= S 2 



"0 

'2 



(5) 



where the reflection amplitude fjK is related to the other amplitudes by tm = —r* M du / d* M . By 
matching inside region I, we find the amplitudes for reflection and transmission through MGS: 
d M = \t\ 2 A 2 ^/(l - A 2>1 A 2tl \r\ 2 ), d_ M = |i| 2 A 2 ,i/_(l - ^2,iA 2 ,i|r| 2 ), r M = r(l - A 2 ,iA 2 ,i)/(l - 
^2,1^2,1 \r\ 2 ), and vm = r*(l — A 2 ,iA 2 , i)/(l — A 2 .\A 2t \ \r\ 2 ). Here we have introduced the amplitude 



Ai.n = Ai(E n , 6) for Andreev reflection at energy E„ 
factors u and v as 



which is defined in terms of the BCS coherence 



Ai(E,t 



Vi(E,6) _E-sgn(E)y/E*-A i (0)* 



Ui(E,9) 



\(0) 



\E\ > |Ai(0)| 
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Fig. 3. Mapping of the process in Fig. [2j on transport through a one-dimensional potential barrier structure in 
energy-space. In (a) the complete structure of the scattering state is shown, collecting the reflections from outside 
regions I-III into the reflection amplitudes r'o and Ty. In (b) the transport through MGS has been collected into a 
scattering matrix. 



Vi (E,0) E-iy/AiW-E* lFl<lAm (a 

where \ui\ 2 + \vi\ 2 — 1 and i = 1, 2 refers to the left and right superconductors respectively. The 
barred amplitudes A are calculated at the angle 6 = tt — 8. 

The amplitudes for reflection from energies outside regions I, II, and III, all the way to ±oo 
in energy (taking into account all possible processes) can be collected into tq and ri respectively. 
The exact numerical values of r§ and are determined by solving the matching equations for the 
complete scattering state by the method described in Ref. Jl(| . Without solving these rather complex 
equations analytically, we can write down the following formal relations between coefficients within 
part II and III respectively: 

bo=r c^ d%=r 2 a%. (7) 

At the injection and exit points, we get four equations connecting part I with the parts II and 
III: 

Oq = JS ae ^ + Aifibg, eg = J8<jh^ + Aifido 

where J = [u\{E) 2 — vi(E) 2 ]/ui(E) is the amplitude for injection into the normal region from the 
left superconductor. The S- functions in Eq.(|^) are included because we do the matching for incoming 
electron-like (a — =>• S ae ^ — 1) and hole-like (er = 8 a h^ = 1) quasiparticles separately. 

The two contributions are summed up in the end in the current formula. The contributions to the 
two-particle current from particles injected from the right superconductor will be proportional to 
V 2 since the trajectories are non-resonant (they do not hit MGS since it only appears at the surface 
of the right superconductor) . We therefore neglect these contributions to the current in the present 
discussion. The structure of the above matching equations is drawn in Fig. ^p. 

Solving Eqs. (§)-(§), we get the following expression for the current density by inserting the 
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obtained coefficients a n , 6„, c n , and d n (n = 0, 2) into Eq. (|4|): 

E 



^r(0,E)K £ r(^^)=2(l-4 )(l-^ 2 )|G 2 o| 2 (l + |r 2 |X 2 )(l + kolXo), (9) 



a- cr—e~',h~' 

where we introduced 

G 20 = ^ ^. (10) 

(1 - A\ 2 r 2 r M ){l - A\ r r M ) - Al Al 2 r r 2 d M dM 

The factors in the current density given in Eq. @ can be interpreted in a suggestive way (see 
also the more general discussion in Ref. (26)). The factor of two appears since we are considering a 
two particle process (transfer of the charge 2e). The factor (1 — A\ ) is the probability to enter the 
normal region at energy Eq and the factor (1 — A\ 2 ) is the probability to exit at energy E 2 . The 
propagator G 2 q is taking us from the injection point up to the exit point. The factor (1 + \r 2 \ 2 A 2 2 ) 
adds the two possible ways of exiting the normal region. First, we may go out directly when we 
come from below (giving the term 1). Second, we may Andreev reflect at E 2 , go further up and 
be reflected back (in region III) and then go out at E 2 (giving the term \r 2 \ 2 A\ 2 ). The last factor 
(1 + ro| 2 ^4i q) adds the contributions from injected electron-like quasiparticles (giving the term 
1) and hole-like quasiparticles (giving the term |ro| 2 A^ ). The reason for getting the extra factor 
|ro| 2 ^4f o f° r the hole-like quasiparticles is that the injection is downwards [into the coefficient cq as 
seen in Eq. (||)] meaning that the particle must be reflected at negative energy, and then be Andreev 
reflected at Eq before going up in energy and out from the normal region at E 2 . 

Eqs. (^) and ([!(]) describe the current density from quasiparticles incoming from the left electrode 
at energy E w —eV, resonantly transmitted through MGS, leaving the normal region into the left 
electrode at energy E 2 w eV. This process is not possible for voltages below eV — A s , because then 
cither E or E 2 is within the gap of the left superconductor. For voltages well above eV — A s we may 
approximate the propagator G 2 o by the bare transmission amplitude du-, since A\$ ss A\ t2 « 0. For 
the particular orientation a = 7r/4, |c?m| 2 is equal to 1/(1 + E 2 /T 2 ), i.e. a Breit-Wigner resonance of 
width r = £>|A(#)|/[2(1 — T>)]. The total contribution to the current, after integration over energy, 
is then I oc T, i.e. / is proportional to 2?, which is of the same order of magnitude as the single 
particle current. This means that the main current (current proportional to T>) has an onset at 
eV = A s , due to the resonant two-particle process. Fig. |J confirms this picture: we see that the 
resonant two-particle process (solid line) gives the main contribution to the total current (dot-dashed 
line) for voltages between A s < eV < Ao. For voltages eV > 2A S the single particle contribution 
is noticeable. We note in passing that the above effect is similar to the effect of a Breit-Wigner 
resonance in the normal region of an s\s junction f28f| . 

Close to the main current onset at eV = A s , the magnitude of the Andreev reflection amplitudes 
|-<4.i,o| an d 1-^1.2 1 are both close to unity. This together with the reflections from ±oo in energy (given 
by ro and r 2 ) give rise to an additional resonance in G 2 q (the denominator gives rise to a singularity 
at the gap edge) . This boundary resonance is quite weak and does not change the order of magnitude 
of the current on its own, since it is connected to the divergency in the superconducting density of 
states which is integrable. At the onset (eV = A s ) the boundary resonance overlaps with the MGS 
resonance. The overlap broadens the bare transmission amplitude which results in a current peak. 
The boundary resonance is due to the reflection from ±oo in energy, so by removing them by hand 
(putting ro = r 2 = 0) we find that the current peak dissappears and is replaced by a smooth onset 
(the dashed line in Fig. ||). 

We note that the N|d w /4 case is easily reached by letting A s — > 0. In this limit the current onset 
will be at eV = and the current peak disappears since we lose the boundary resonances (ro and 
r 2 -> when A s -> 0). The onset at eV = result in a ZBCP &p|. 
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Fig. 4. Comparing the current contribution from the two-particle process (solid line) with the complete IV-curve 
(dot-dashed line), we see that the peak is due to this particular process. For comparison, we include the current 
(dashed line) from the bare MGS resonance (letting ro = V2 = 0). Taking into account the reflections ro and T2 gives 
rise to boundary resonances, which at the voltage eV = A s = 0.2Ao overlap with (and therefore broadens) the MGS 
resonance and produce the current peak. 



4. Zero-bias conductance peak in the d a \d- a junction 

In recent experiments on grain boundary junctions of HTS a ZBCP has been found [^o[ Ell. 
Within our model a ZBCP can only be found in rather transparent d a \d- a junctions as previously 
reported []l0| . Generally, the current contribution from a process of order n is (without resonances) 
proportional to the transparency T> of the junction to the power of n |24|, [2f| . In the limit of small 
bias voltage, the current contributions are from high-order processes, meaning that the current is 
extremely small without any resonances. This picture is changed by the MGS resonance. 

In Fig. H we have drawn the map of the scattering state contributing to the current at small bias 
voltage. The injected quasiparticles undergo multiple Andreev reflections and exit the normal region 
above the gap. At some particular voltage, the lowest order process contributing to the current, will 
contain 2n Andreev reflections and MGS will be reached, after n Andreev reflections, at the sideband 
energy E n . The path from the injection point up to MGS consists of n passings through the barrier. 
In Fig. [s]b we have collected this path into an effective first-order path with a barrier of height V n . 
The same is done for the path from MGS up to the exit point. The map in Fig. |5|b is analagous to 
the map in Fig. ||a, meaning that we should expect a Breit-Wigner resonance of width Te// oc D n , 
broadened by the boundary resonances. In Fig. ||c the middle region (corresponding to region I 
in the discussion of the s\d^/^ junction) is collected into a scattering matrix. We can from this 
picture immediately write down a formal expression for the current, equivalent to Eqs. (||) and (fl0|), 
in terms of the effective reflection and transmission amplitudes through MGS {reff, ^eff, ^e//j 
and deff), and the reflection amplitudes from ±oo in energy (ro and T2 n )- Since the process under 
consideration is of order 2n, the expression for the current will contain the factor 2n. Due to the MGS 
resonance, the current density will therefore contain peaks of height 2n and width Te// , which after 
the integration over energy give rise to an enhanced current at small voltage resulting in a ZBCP. 
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Fig. 5. Mapping of the 2n particle process involving MGS contributing to the current in the high transparency 
d a \d— a junction at small voltage. The n barriers on each side of MGS shown in (a) are in (b) collected into effective 
barriers with transparency X> n . The resonant transport through MGS can then be described by an effective scattering 
matrix as drawn in (c). The bare transmission resonance (described by rfe//) is broadened by boundary resonances 
(due to the reflections rrj and r^n at the boundaries) producing an enhanced current at small voltage, i.e a ZBCP. 



Since the peak width scales with V n (neglecting broadening effects due to boundary resonances), 
the current contribution will be considerable only in junctions with high transparency JTo| ]. 



5. Summary 

We have discussed the effects of the midgap state on current transport through voltage-biased 
Josephson junctions of d-wetve superconductors and demonstrated that in connection with multiple 
Andreev reflection the MGS acts as a transmission resonance in energy space. 

Depending on the orientation of the d-wave superconductors and the transparency of the junction, 
resonant transport through MGS influences the current in different ways. In low transparency s\d n /4 
junction the MGS resonantly enhances the two-particle current, which gives the main contribution 
to the current for voltages between A s < eV < Ao. At the onset an overlap of boundary resonances, 
connected to the peaks in the density of states, and the MGS resonance results in a current peak. 
In the same way, the MGS resonantly enhances the two-particle current in d\d- junctions which, 
depending on the orientations of the superconductors, results in a current peak at finite voltage Q . 
Since the gap has nodes in both electrodes there is no clear onset of the resonant current. In addition, 
a smooth background contribution from the single-particle current at all voltages makes the effect 
less pronounced compared to the sjc^^ junction. The relation between our results and experiments 
was discussed in Ref. JlO[ . 
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In high transparency d a \d- a junctions there are resonant 2n-particle processes involving the MGS, 
corresponding to tunneling through an effective symmetric double barrier structure in energy space. 
These symmetric processes give a huge contribution to the current, resulting in finite current at 
low voltages, as in ballistic s|s-juntions. This may explain the zero-bias conductance peaks seen in 
experiments on grain boundary junctions of HTS superconductors 20, 2l| . 
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